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Abstract The energy levels of the double-well potential receive, be- 
yond perturbation theory, contributions which are non-analytic in the 
coupling strength; these are related to instanton effects. For example, 
the separation between the energies of odd- and even-parity states is 
given at leading order by the one-instanton contribution. However 
to determine the energies more accurately multi-instanton configu- 
rations have also to be taken into account. We investigate here the 
two-instanton contributions. First we calculate analytically higher- 
order corrections to multi-instanton effects. We then verify that the 
difference betweeen numerically determined energy eigenvalues, and 
the generalized Borel sum of the perturbation series can be described 
to very high accuracy by two-instanton contributions. We also cal- 
culate higher-order corrections to the leading factorial growth of the 
perturbative coefficients and show that these are consistent with an- 
alytic results for the two-instanton effect and with exact data for the 
first 200 perturbative coefficients. 
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In theories with degenerate minima, the energy eigenvalues of the states at nonvanishing couphng 5 7^ 
cannot in principle be obtained by analytic continuation from the unperturbed situation at vanishing 
coupling g = because a potential with degenerate minima introduces a degeneracy in the spectrum: for 
any one unperturbed state, two states emerge when the perturbation is switched on. These two states 
are seperated by an energy shift which is nonperturbative and nonanalytic in the coupling, i.e. vanishing 
to any order in perturbation theory. Therefore, the two states are described by the same perturbation 
series and yet differ in their energy by instanton contributions. Specifically, we consider the case of the 
double-well potential with the hamiltonian 



H--l§^ + lvii) vi,)^l,m-,?. (1) 



It has been conjectured that an asymptotic expansion for the energy eigenvalue can be obtained by 
finding a solution to the equation 



1/1 , A / 2^^(^ 



— ^^^-D{E,g)ji^--j exp[-^(£;,5)/2]=±i, (2) 

which can be understood as a modified Bohr-Sommerfeld quantization condition. The plus and minus 
signs apply to even- and odd-parity states, respectively. The conjecture (|^), whose validity has been proven 
in has found a natural explanation in the framework of Ecalle's theory of resurgent functions 
The functions D{E,g), A{E,g) constitute power series in both variables. The function D{E,g) describes 
the perturbative expansion; its evaluation is discussed in The first terms read 



D{E,g)=E + g + \^ +9^ {ib E^ +'^-^ E^ + 0{g^) . (3) 

The ground and the first excited state are both described by the same perturbation series which can 
be found by inverting the equation D{E,g) = 1/2. When the energy is expressed in terms of the naive 
perturbation series in g, the function D{E,g) then vanishes in any order of perturbation theory, i.e. in 
all orders in g. For a general state, D{E,g) = -I- 1/2 where N is the quantum number of the unper- 
turbed state which is a harmonic oscillator eigenstate. The function A{E, g) essentially describes instanton 
contributions |4j; its first terms read 

A{E,g)^j-+g ^l7E'+^^^+g' ^227 E' + ^-^ E^ + 0{g') . (4) 

A solution to the equation (|2|) can be found by systematically expanding the energy eigenvalue E(g) in 
powers of g and in the two quantities 



9/ vT^g 

Terms of order ^(g)" belong to the n-instanton contribution. The energy eigenvalue for nonvanishing 
perturbation g ^ can be described by two quantum numbers: the unperturbed quantum number N and 
the positive or negative parity of the state. We have (the upper index denotes the instanton order) 

00 

i?^,±(5)-E<±(5) (5) 

n=0 

where the perturbation series (zero-instanton contribution) is given as 

00 

K=0 
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where the right-hand side is parity independent. For n > 0, the instanton contribution reads 



N,±i9) = [ - 



Nn 



nkl y ' 
fc=0 1=0 



(7) 



The lower indices n, k and I of the e coefficients denote the instanton order, the power of the logarithm 
and the power of g, respectively. The results relevant for the current investigation read, 
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where 7 = 0.57221 ... is Euler's constant. Odd-instanton contributions have opposite sign for opposite- 
parity states and are responsible, in particular, for the energy difference of the ground state with quantum 
numbers (0,-1-) and the first excited state with quantum numbers (0,—). The dominant contribution to 
the seperation of the two lowest energy levels is given by the one-instanton contribution: 
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g' + 0{g') +Oiagr). 



(9) 



By contrast, even-instanton contributions have like sign for opposite-parity states and are responsible, 
in particular, for the displacement of the mean value (1/2) [Eo,-{g) + Eo,^{g)] from the value of the 

generalized Borel sum of the perturbation series B ^X]^=o ^o^k 5^) i^'^^ evaluation of the generalized 
Borel sum of a nonalternating divergent series see for example Sec. VI and Table III of ||]). The dominant 
contribution to the displacement comes from the two-instanton effect, and we have 



,(5)+i?o.+ (5)]-Re<^S 
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The function 



1 



A(5)=4. 



[Eo,^{g) + Eo,+ ig)] - Re {b {Zk=, E^% g^) } 

[EoAg)-EoMg)f in(^) 



(11) 



relates the multi-instanton contributions to the energy eigenvalues, which can be evaluated numerically, 
and to the (generalized) Borel sum of the perturbation series which is evaluated by analytic continuation 
of the integration path into the complex plane (see ||l^). The calculation of A{g) at small coupling is 
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problematic because of severe numerical cancellations. From the equations 
the following asymptotics for A (5), 



oh and rtll), we obtain 
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If we additionally perform an expansion in inverse powers of ln(2/g) and keep only the first few terms in 
{l/ln(2/g)} in each term in the ^-expansion, the result reads 



A{g) ~ 1 + 3g - 
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The higher-order corrections, which are only logarithmically suppressed with respect to the leading terms 
1 + 3(7, change the numerical values quite significantly, even at small coupling. In Table we present 
numerical results for the function lS.{g) at small coupling; these are in agreement with the first few 
asymptotic terms listed in equation ( p^ up to numerical accuracy. Of course, for strong coupling, sig- 
nificant deviations from the leading asymptotics must be expected due to higher-order effects; these are 
indeed observed. For example, ai g = 0.1 the numerically determined value reads A(O.l) = 0.87684(1) 
whereas the first asymptotic terms given in equation (O) sum up to a numerical value of 0.86029. 



Table 1: Comparison of numerical values for the function A((7) defined in equation ( |11| ) in the 
region of small coupling to values obtained by calculating the first few terms in its asymptotic 
expansion given in (|l2|). 



coupling g 


0.005 


0.006 


0.007 


0.008 


0.009 


0.010 


A ((7) num. 


1.0063(5) 


1.0075(5) 


1.00832(5) 


1.00919(5) 


1.00998(5) 


1.01078(5) 


A(g) asymp. 


1.00640 


1.00739 


1.00832 


1.00919 


1.01001 


1.01078 



The higher-order corrections to the two-instanton effect are related to the corrections to the leading 
factorial growth of the perturbative coefficients. This can be seen by expressing that the imaginary part 
of the perturbation series, when continued analytically from negative to positive coupling, has to cancel 
with the imaginary part of the two-instanton contribution which is generated by the logarithms \n{—2/g). 
The corrections of order g ln(— 2/5) and g^ ln(— 2/g) yield the \/K~ and l/X^-corrections to the leading 
factorial growth of the perturbative coefficients. From the results for ^''2if'^ {j = 0, 1, 2) given in equation 
(H), we obtain 
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The analytic results should be checked against explicit values of the perturbative coefficients. We have 
determined the first 200 perturbative coefficents e'^]^ (K = 0, . . . , 200) of the perturbation in the form of 
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rational numbers, i.e. to formally infinite numerical accuracy. This allows to verify the 1/K- and 1/K^- 
corrections to the leading factorial growth in equation ( |l^ to high accuracy, for example by employing 
Richardson extrapolation Using the 160th through the 200th perturbation coefhcient as input data 
for the Richardson algorithm, the coefficients of the leading, of the l/_fi'-subleading and of the 



suppressed corrections are found to be consistent with the analytic results given in equation ( 14 ) up to a 
relative numerical accuracy of 10^^^, 10^^'^ and 10"^", respectively. For completeness, we give here the 
numerical values of the 198th through the 200th perturbative coefficients, to 30 decimals. These read: 

£:(°)gg = 5.50117 76962 885879352775694 38632 x lO''*'" , 
E^°lgg = 3.28445 39841 65780 00616 21912 32835 x 10*^^ , 

■^o°2oo = 1-97082 14193 09543 76979 53006 07410 X 10^™. (15) 

Values for all 200 coefficients will be presented elsewhere. 

It is an interesting consequence of the expansion that the energy difference (i?o,- — -E'0,+)j at small 
coupling, is described to high accuracy by the one-instanton contribution (n = 1 in equation (^). For 
g — 0.001, we obtain to 180 decimals, 

£;o,+(0.001) = 0.49899 54548 6210917168 91308 3948192163 68209 47240 
20809 66532 93278 69722 01391 15135 28505 38294 45798 
45759 95999 06739 55175 84722 67802 81306 96906 01325 
25943 77289 94365 88255 24440 17437 12789 27978 99793 , (16) 



whereas 

i;o,_ (0.001) = 0.49899 54548 6210917168 91308 3948192163 6820947240 
20809 66532 93278 69722 01391 29839 92959 55803 70812 
27749 92448 48259 36743 64757 68328 84835 3551 1 34663 
06309 82331 51885 23308 08622 84780 52722 10103 67282 . (17) 

Decimals which differ in the two energy levels are underlined. The results have been obtained by lattice 
extrapolation using a modified Richardson algorithm which is constructed according to ideas outlined 
in [Q. Calculations were performed on IBM RISC/6000 workstations while making extensive use of 
multiprecision libraries [p"3|-p^. We define Vnig) as the Mth partial sum of the one-instanton contribution 

<1(5) -<!(<?), 

M 

VM{9)^2^{g)Y,e^°,7^g^. (18) 

Using exact rational expressions for the coefficients Cioj (j < 141), we obtain 
Pi4o(0.001) X 10^1 = 

1.47046 44541 75092 50138 19899 64494 15198 15678 00350 05260 35283 , 
86053 33378 03660 50415 75193 50528 41826 73433 99328 21246 74888 (19) 

Pi4i(0.001) X 10^1 = 

1.47046 44541 75092 50138 19899 64494 15198 15678 00350 05260 35283 , 
86053 33378 03660 50415 75193 50528 41826 73433 99328 21246 74887 . (20) 
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These values are in excellent agreement with the numerically determined energy difference (see the results 
presented above in equations ( [f6| ) and (|f^)) 

[So,-(O.OOf) - £;o,+(0.00f)] X fO^i = 

1.47046 44541 75092 50138 19899 64494 15198 15678 00350 05260 35283 , 
86053 33378 03660 50415 75193 50528 41826 73433 99328 21246 74887 . (21) 

The first 70 decimals in equations ( |l6|) and (^) are the same because the one-instanton contribution is 
of the order of 1.4 x 10^^^. The accuracy to which the one-instanton contribution describes the energy 
difference i^o,- (0.001) — i?o,+ (0.001) is limited by the three-instanton effect which for g = 0.001 is of 
the order of 8 x 10^^^^. Note that the two-instanton effect (which for g = 0.001 is of the order of 
4 X lO^^^'^) does not limit the accuracy to which the one-instanton contribution describes the energy 
difference because it has the same sign and equal magnitude for opposite-parity states with the same 
unperturbed quantum number N. 

We have demonstrated that the behavior of the characteristic function A{g) defined in equation ( ITTI ) 
at small coupling is consistent with higher-order corrections to the one- and two-instanton contributions, 
specifically with the instanton expansion of the energy levels governed by the equations (^ and (^ , with 
the assumption that the instanton contributions given by equation (^ should be Borel summed, with the 
explicit results for the higher-order coefficients listed in (||) and the analytically derived asymptotics for 
the function A{g) given in equation (|l2|). The corrections of relative order to the leading factorial 

growth of the perturbative coefficients - see equation (|l^) - are consistent with the analytically evaluated 
(7™ ln(— 2/(7)-corrections to the two-instanton effect and with the explicit values for the first 200 terms 
in the perturbation series (|^). The nonpcrturbative energy difference _Bo,-((?) — i?o,+ (5) at small coupling 
g is described, to high accuracy, by the one-instanton contribution only. 
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